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Introduction 



Matrix models [3', play a very special role in modern theoretical physics. They appear regularly 
and prove useful in analysis of various simplified models of concrete physical phenomena, but their 
real significance is that they somehow capture and reflect the very basic properties of string theory 

- and can serve to represent the universal classes of quantum field theory models. From the very 
beginning matrix models were introduced to describe some very general features (eigenvalue repul- 
sion) of statistical distributions [5]. It was, perhaps, the first recognition of the role of group theory 

- the underlying theory behind matrix models - in explaining the fundamental properties of quan- 
tum/statistical behavior. Much later this led to discovery that integrability is the basic property of all 
functional integrals, considered as functionals on the moduli space of theories [6], and matrix models 
played a central role [7J in the formulation of the fundamental relation 



partition function = r-function (1) 



between the two central concepts of modern theory, already with a variety of applications in different 
fields, from gauge theories [8] to Hurwitz theory [9] and with still many more to come. An immediate 
implication of ([T|) is that quantitative approach - a possibility to calculate something - in string 
theory (= a theory of families of quantum mechanical models) requires extension of the standard 
set of special functions to a broader set of r-functions \10\ [TT] - a far-going generalization of both 
hypergeometric and elliptic families. A highly non-trivial step here was introduction of " infinite-genus" 
T-functions, satisfying the string equations [12] and Virasoro/VF-constraints |13j-|21j. and it was once 
again inspired by the study of matrix models. Unfortunately, even the simplest of these r-functions, 
associated with Hermitian |3j and Kontsevich |22 1 l23 1[Mll25j matrix models, are not yet systematically 
studied/tabulated and still can not be included into the special-functions textbooks ~ see [26] for the 
first attempts in this direction. It is very important to realize that the world of such r-functions is 
cognizable, and is, perhaps, actually finitely-generated: many (all?) matrix-model r-functions are 
actually expressed by group-theoretical methods through a few basic ones. This decomposition results 
from description of genus expansion of matrix-model partition functions in terms of auxiliary "spectral" 
Riemann surfaces [27 ^1261 [28] . and explicitly relates them to representation theory of Krichever-Novikov 
type deformation [29] of Kac-Moody algebras. One of spectacular byproducts of this development is 
the possibility to build a " string- field-theory- like" diagram technique [28j for the model of entire 
string theory, provided by "M-theory of matrix models" f2]. As shown in [1], the main basic block 
(constituent) of the matrix model partition functions in this approach is the ordinary Kontsevich 
r-function Zk of p2j . 

However, already in pLj a first counter-example was found to this (over?)-optimistic conjecture: 
partition function of the complex matrix model [30[ [3T] is made not only from Zk , but also from some 
other ingredient, denoted Zk in s.8 of [I]. The purpose of the present paper is to identify this Zk with 
a very important and well-known partition function: that of Brezin-Gross-Witten model (BGWM) 
[3211M112D]: 

Zk = Zbgw (2) 



By definition, BGWM describes correlators of unitary matrices with a non-linear Haar measure. Uni- 
tary correlators play a crucially important role in description of gluons in lattice gauge-theory models 
[34] . however, unitary matrix models are more complicated than Hermitian ones, they are in inter- 
mediate position between eigenvalue and non-eigenvalue models and remain under-investigated, see 
p5]-f38] for some crucial references. A modern matrix-model-theory approach to BGWM and its 
embedding into the set of generalized Kontsevich models (GKM) [23] was outlined in [20], but has 
not been developed any further since then. Hopefully reappearance of this model in the context of 
matrix-model M-theory will help to attract new attention to this unjustly- abandoned subject. 

We begin in sUfrom reminding the definition and the main properties of the four partition functions 
which participate in our story: Znit), Zc{t), Zk{t) and ZsGwij)- They were originally introduced 
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as matrix integrals over Hermitian {Zh and Zk), complex {Zc) and unitary (Zbgw) matrices, with 
the time-variables identified either with the coupling constants {t^ in Zh and Zq) or with the Miwa 
transform of the background field (r^ in Zk and Zbgw)- 

As functions of their parameters - the time- variables t or r - these integrals satisfy Ward identities 
(or Picard- Fucks equations) [T^, which have the form of the Virasoro constraints. Namely, 

— = Lk^2ZH, -TTT^ = Lk-iZc (3) 
dtk otk 

with k > 1, with = N/gZn, = N/gZc, where A'" is the size of the matrix in the original 
integral representation and with " discrete- Virasoro" operators |15j 



d ™ (?2 



a+6=m 

and 

dZx ^ dZsGW ^ 

7^ = L-k-2ZK, = Lk-lZBGW (5) 

OTk-1 OTk-1 

with k >1 and with " continuous- Virasoro" operators [1'6\ 114] 

°° / l\ d d"^ 1 

V 2 J dTm+k dTadn I6g^ 16 ^'^) 

a+b=m—l 

Now one can switch to D-module approach and define the four partition functions as solutions to the 
four systems of linear differential equations ([3]) and ([5]), and original integral formulas are just integral 
representations for the solutions. 

Not surprisingly, such representations are not unique, and one can instead represent the same 
solutions in a very different integral form: of Kontsevich-Penner integrals over nxn Hermitian matrices 
with a peculiar Penner term A^tr log^ in the action. This puts all the four models in the unifying 
context of GKM theory [53]. Direct relation between the two integral representations is provided by 
a version of Faddeev-Popov trick from [21] . 

All four matrix integrals can be expressed in the form of determinants of the other matrices, which 
have ordinary single integrals as their elements. These determinant representations are very important, 
because they are typical for the tau-functions of integrable hierarchies [401 [24] - the generalized 
characters of Lie algebras [39|lll]. In other words, partition functions of the matrix models are always 
the tau-functions [7J. Moreover, this is a general property of all partition functions - the generating 
functions of all correlation functions in any quantum theory, - this is a consequence of the freedom to 
change integration variables (fields) in the functional integral [6j. For tau-functions the whole sets of 
Virasoro constraints are actually fixed by their lowest components L_i or Lq, which therefore has its 
own name: the string equation |12] . 

Integrability means that partition function satisfies a bilinear Hirota equation [10] of the form 

This equation has its origin in decomposition rule Rx R' = J2i ^or representations of Lie algebras 
and this is an equation for the characters of the algebra [69\ [TT] . For loop algebras the characters can 
be rather non-trivial, they can be actually labeled by some auxiliary (spectral) Riemann surfaces or, 
better, by the points of an infinite-dimensional Grassmannian (the universal moduli space of [41]) - 
what means that the spectral surface can actually have an infinite genus (and this is typically the case 
for the matrix- model partition functions). 
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For irreducible representations of finite-dimensional simple Lie algebras, the characters are given 
by two determinant Weyl formulas: in the case of SL(N) and representation labeled by partition 
m : rrii > m2 > • • • > rrifyf > (the Young diagram) 

Xm{t) = detPi+m^^i{t), (8) 
where Pi{t) are the Shur polynomials, exp (j^k'^kx'^^ = J2m -^mit) , and, after the Miwa transform 

In fact, the two Weyl formulas are mirrored in two possible representations of matrix models 
partition function: as we shall demonstrate, the Hermitian and complex matrix models, besides the 
standard determinant representation of type ([8]), have a Kontsevich-Penner representation of type ([9]). 

Virasoro constraints can be used as recursion relations to provide the logarithm of the partition 
functions logZ in the form of the formal series in non-negative powers of t- variables and g^. Such 
formal series are unambiguously defined by the systems ([3]) and ([5]), since A; > 1 in all the four cases. 
The generating functions (" multiresolvents" ) p^^'^\z), defined as 

oo 

...,zg)= V(zi) . . . V(z,) log = ^ g^P-^p<P\'^\zu ...,z,) (10) 

p=0 

possess an important property: they are poly-differentials on auxiliary spectral Riemann surfaces 
(complex curves) [26l, 128], which for the four matrix models in question are all double-coverings of 
the Riemann sphere with only two ramification points (and thus Riemann spheres themselves). The 
spectral curve representation arises only for the special choice of generating functions: they should be 
resolvents, i.e. no /c-dependent coefficients are allowed in (jlOp . There are other interesting choices of 
coefficients, when alternative generating functions possess other interesting properties, see for example 

m- 

In sin] we proceed to decomposition formulas. Virasoro constraints can be also considered as 
quadratic differentials on the spectral surfaces, expanded near particular points. It turns out that 
"discrete" operators arise in expansion near non-singular points, while "continuous" operators ([6]) 
- in those near ramification points of degree 2. This means that a globally-defined Virasoro quadratic 
differential can be decomposed in both bases and this idea finally leads to decomposition formulas [1]. 

The basic one is 

ZH = UKKiZK^ZK), (11) 

it expresses Zfjit) for Hermitian matrix model through Zc{t) for Kontsevich model. It is explained in 
full detail in sl3l Ingredients of the construction are: explicit parametrization of the spectral curve and 
of the singular differentials in the vicinities of particular points on it, explicit formula for the global 
U{\) current and the Virasoro differential, its projection onto "canonical" quadratic differentials in the 
vicinities of the particular points and Bogoliubov transform of the time-variables with the help of the 
conjugation [/-operator, and, finally, projection from generic Laurent series for the Virasoro operator 
to a Taylor series, provided by peculiar projection operator V, which picks up a triangular subalgebra 
from entire Virasoro (Krichever-Novikov) algebra. Generators of this triangular subalgebra can be 
imposed as constraints on partition function and form a consistent and resolvable set of constraints. 

Prom the point of view of D-module approach the difference between Zh and Zq in ([3]) looks 
minor: both are annihilated by the same discrete- Virasoro operators L„; only n > — 1 for Zh, but 
n > for Zc- The second difference is that the shift of time- variables, which generates the l.h.s. 
in ([3]), is also different: t2 is shifted in the case of Zh, but ti is shifted in the case of Zc ~ this is 
important to explain why both sets of equations ([3]) have unambiguous formal-series solution, despite 
the set of constraints contains one less equation in the case of Zc- However, this second difference 
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is not essential for comparison of ([3]) and ([5]). Indeed, the relation between Zk and Zk = Zbgw 
in ([5]) is exactly the same: both partition functions are annihilated by the same continuous- Virasoro 
operators Cn, but n > — 1 for Zk, while n > for Zbgw- This time shifted are ri in Zk and tq in 
ZbgWj what guarantees that the formal-series solutions are unambiguously defined in both cases. All 
this implies that in the spectral-surface formalism the difference between Zfj and Zc is concentrated 
in the choice of projection operator V at the last stage. Since projection operator can be realized as a 
contour integral, its modification can actually be shifted from Zh (where it transformed Zjj into Zc) 
to one of the two Zk at the r.h.s. of (jlip and convert it into Zk = Zbgw ■ In other words, instead of 
PT]) we obtain 

Zc = UKk^ZK®ZK) i.e. Zc = tjj^i^{ZK®ZBGw) (12) 

This is the main result of the present paper and it is discussed in full detail in slU This formula was 
supported in [T] by explicit comparison of the first terms of expansions of partition functions at both 
sides of (fT2]) . 

Some concluding remarks are given in slH 
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Part I 

The four matrix models 



1 Hermitian matrix model Znit) 
1.1 Original integral representation 

Partition function of Hermitian matrix model [3] is defined by the integral over hermitian N x N 
matrix 

Znit) = I exp ( -T^Tr H'' + - V t^Tr hA dH (I.l) 



where the measure dH = JJj^ j^^dHij. This is nothing but the generating function of all GL{N)- 
invariant Gaussian correlators of Hermitian matrix H. 

1.2 Eigenvalue representation 

The Hermitian matrix $ can be diagonalized by a unitary transformation, H = UDU\ where D = 
diag(ifi, . . . , Hj\f) is the diagonal matrix made from the eigenvalues of H. The norm of H decomposes 
as 

2 N N 

Tr {6Hf = Tv {6Df + Tt {ju^6U,D]) =Y,{5Hif + Y,{Hi - Hjf{U^ 5U)l (1.2) 

1=1 i<j 

SO that the measure 

N N 

dH = [dU] WiHi - Hjf W dH, (1.3) 

where [dU] is the non-linear (Haar) measure for unitary matrices. Since the action in does not 
depend on U, the integral Vat = /[dC/] fully decouples, and (jl.ip turns into an A^-fold integral over 
eigenvalues H, with the peculiar square of the Van-der-Monde determinant A(/7) = Y\f--^j{Hi — Hj) 
in the measure: 

ZH{t) = VM j\{dH,^\H)e^jyl-^Y.H^ + -Y.tkHi\ (I-^) 

i \ ^ i ^ i,k J 

1.3 Virasoro constraints 

Integral p.ip is invariant under any change of integration matrix- variable H. In particular, Znit) does 
not change if one substitutes H ^ H + eH^^^ with any matrix-valued parameter e and any integer 
n > — 1. This invariance implies that [16] 



^^^^^^ - KZnit), n>-l (1.5) 



dtn+2 

with the operator Ln defined in dH) and 

oto g 

The l.h.s. in ()I.5P is produced by the shift of the t2 variable ^2 ^2 ~ 1/2 in the initial formula p.ip . 

Together with 0.6^ the system O.SP provides a set of recurrent relations which allows one to 
unambiguously construct Zuit) term-by-term as a formal series in non-negative powers of t- variables. 
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1.4 Determinant representations and integrability 

As we shall explain now, the properly normalized matrix integral Zjq = ^ ^^ Zfiit) is a r-function of 

the Toda-chain integrable hierarchy [15^ \7\ . In this paragraph we rescale the time variables to cancel 
the coefficient g in front of them in order to have the standard definition of integrable hierarchies. 

One of the technical ways to deal with integrals of form ()I.4|) was proposed in [3]. The authors 
introduced a system of orthogonal polynomials with the orthogonality condition 

J P,{H)Pj{H)e-''^''^dH = (^,,e'^"W, (1.7) 

where e'^'^*-' are norms defined by integral ()I.7P and the normalizing condition for the polynomials is 

PiiH) = ^-iijW, -in = 1, (18) 

i.e., the coefficient of the leading term is put equal to unity. 
Using polynomials (\l.7\i . (jl.Sp . we can rewrite (jL4]) as 

1 

Zm = (Niy' J I[dHidetPk-iiH,)detPi.iiH^)exp{--J2Hf + J2t^H^} = J] e^'^'l (1.9) 

i i i,k i=0 



In order to get the determinant representation, we rewrite orthogonality condition ()I.7p in a "ma- 
trix" form. That is, we introduce the matrix T with matrix elements ^mn 

defined by (|I.8P , the so-called 

moment matrix C with the matrix elements 

dj = J dHW+^-'^e-^^^\ (1. 10) 

and the diagonal matrix J with diagonal elements e'''" . Then, (jLT]) can be written as a matrix relatioE0 

rcT'^ = J (i.ii) 

where V'^ is the transposed matrix. Evaluating the determinant of the both sides of this relation and 
using (|I.9p . we obtain 

Zn = det dj. (1.12) 

Ny.N ■' ^ ' 

The moment matrix satisfies a number of relations, which follow directly from its explicit form p.lOp . 

^ - dkC^t) = ^ a'=a(t), (1.13) 

a, = c,+,- (1.14) 

and 

= a^-2cn = (1.15) 

Finally, the partition function of the one-matrix model is 

Zn = deid'+^-^C, (1.16) 

which results in the Toda chain [21 SOI [3 |15] . Note that conditions (lLT3ll and (ILT5|) are satisfied 
for the whole hierarchy of the two-dimensional Toda lattice and for the KP hierarchy, while p.l4p is 
specific for the Toda chain. 



^This relation is nothing but the Riemann-Hilbert problem also known as the factorization problem, see [431 1151 144] 
for the details. 
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1.5 Kontsevich-Penner representation 

The system p.Sp is solved by another integral, very different from p.ip |45| : 

Zj/(t) =7r^^5^^i-^"e^^^ ^ / exp --tr /i^ + iVtr log/i-itr /iL d/i (1.17) 

Jnxn V 2 / 

Integral is now over n x n Hermitian matrix h, dh = na,6=i dhab and depends on additional n x n 
matrix (background field) L. To emphasize the difference between A'^ and n we use small letters for h 
and tr instead of H and Tr in O.lh . If expanded around a saddle-point h = L this integral is a formal 
series in positive powers of variables 

= -|tr L-^ k > 1 (1.18) 

ancH 

to = <7trlogL (1.20) 

This is a model from the GKM family and peculiar logarithmic term in the action is often named 
Penner term [36], so that (|I.17p is known as Gaussian Kontsevich-Penner model. 

1.5.1 Proof I: Ward identities 

In order to check that O.lTh satisfies ([3]) one begins with the Ward identity for this integral [23], 
associated to the shift h ^ h + e oi the integration variable by a small arbitrary matrix e: 

which gives 

It remains to substitute a function Znit) with fs expressed through L by the Miwa transform p.isp . 
Then (OB becomes see [IB |2ll [7] for technical details. 

1.5.2 Proof II: Orthogonal polynomials 

The two integral representations p.ip and p.l7p can be related directly, without a reference to Virasoro 
constraints p.Sp . One procedure, making use of orthogonal polynomials (Hermite polynomials in this 
particular case) is described in details in [24^ [7]. One can rewrite p.OP with time variables substituted 

^It is possible to introduce dependence on to implicitly, namely to define the partition function as follows: 

Z«W=^^<7^^i-^"e^ + ^^'"^'-^*^'°^^ / e^p ( -hr + Ntr log h-^tr hi) dh (1.19) 
this slightly change all calculations but leave the partition function unchanged. 
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by p.isp . (|L20p in terms of orthogonal polynomials in the following way: 

zn = m-' I \{dH,^\H)e^v ( - E^(^^) + - E *.^n 

= I n dH,^\H) exp ( - E V{Hi)\ - Hi) = 

i \ i / i,a 



A(L) 



{N\r^ J n dH, exp E A{H) 
= {N\)-^A-\L) J YldHiew (^-E^(^i)) X 



(1.23) 



X det Pi^iiHA det 

AfxAf (Af+n)x(7V+n) 



Pi-l{La) '■ PN+b-l{La) 



where i,j = 1,...,N, a,b = l,...,n and the orthogonal polynomials Pk{H) are orthogonal with the 
measure exjp{—V{H)). In the case under consideration, it is equal to exp (^—^H^^. Now calculating 
the determinants and using the orthogonality condition p.7p . one arrives at 



Zn = A-\L) det pN+a-iiLb) n e 



•fids) 



n 



det(ab) cl)i^\Lb) 
A{L) 



det(afe) (^i^^(Lfc) 
A(L) 



(1.24) 



with 



<Ai^)(L)=PAr+„_i(L) 



(1.25) 



Now let us see that integral (|I.17p can be also transformed to this form. To this end, we need the 
Itzykson-Zuber formula, [U] 



U(n) 



"A(x)A(y) 



(1.26) 



where integral runs over unitary nxn matrices with the Haar measure dU, and Oj, bj are eigenvalues of 
Hermitian matrices A and B. Now, using p.Sp . we can perform integration over the angular variables 
and rewrite ()I.17p as 



f tr l2 



dhexp [ ~^tr /i^ + A^tr log h — itr hL 



s» / n rf^^ ^ (- f E + ^ E i°g - ^ E ^^^^) 

i ^ ' \ i i i / 



(1.27) 



n '^^^ exp f - ^ E + ^ E log - ^ E ^^^^ 



det <&i(Lj) 
A(L) 



where 



^>i(L) = 629-^" / dxx* ^exp ( -^x^ + TVlogx - ixL 



5 2 



(1.28) 
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It remains to note that the orthogonal polynomials with the measure exp (^—■^H'^^ are the Hermit 
polynomials which have the integral representation 

Pkix) = ^ / dyl^ + iy] e-y^' = / dy/e"^""^ (1.29) 

\/2ttJ-oo \V5 / v2vr J-oo 

This finally reduces (Oijl to (fLTT]) . 

1.5.3 Proof III: Faddeev-Popov trick 

Another way to connect the two matrix integrals, suggested recently in [21j is by using the Faddeev- 
Popov trick. In order not to make calculations with Grassmann variables, we choose the opposite sign 
in (fL18]) . (fL20l) 

4 = |tr L-^ k > 1 (1.30) 

rv 

to = -gtilogL (1.31) 

Then, p.l7|) should be substituted with (the results for the two choices of sign can be also related by 
continuation) 



e-i*'^-^"/ exp(--tr /i^-iVtr log/i + tr (1.32) 

Jnxn V 2 J 

If Miwa transform (|I.30P is made in the original integral (|I.ip . it becomes [7] 

m = (detL)-^ [ , '^J'Z'^^r 1^ = / / /e"^^' (1.33) 
^ — ^ ^ ' JNxNdet{I(g)I-H(g)L^^) J J J 

where Faddeev-Popov trick is applied to substitute the determinant in the denominator by an integral 
over auxiliary rectangular N x n complex matrix fields B and C = . Here (PB = dBdC = 
niliOLi^^'^^aand 

B{H ®I -I® L)C = BaiHijCja - BiaLabCbi = tr B^HB - Tr BLB^ (1.34) 
Taking the Gaussian integral over H we finally obtain: 

dLH) = 5^ / / ef Tr BBtBBt+Tr BLBt^2^ (135) 

At the same time the Kontsevich-Penner integral (jL32]) is equal to 



(jLMl) = g-+- / — — -jr = g 



(1.36) 



det (L + ghr -J det{I gh + I L) 

= 3^+4 1 1 le-ft^ ^'+^i'^>^^+''^^)^'dhdBdC 

where the fields B, C are exactly the same as in p.33p . while 

B{I®gh + I® L)C = gBiahabCbi + BiaLabCbi = gTr BhB^ + Tr BLB^ (1.37) 
Again we can take the Gaussian integral over h and obtain: 

dLMD = 1 lefTr bb^bbUTt blb^ ^2-^ (j_38) 

i.e. exactly the same expression as at the r.h.s. of p.35p . Thus we conclude that 

(HH) = (1-39) 
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the two integral representation for Z}{{t) coincide. 

Inverting the argument, the two matrix-integral representations p.ip and ()I.32p for Znit) are 
associated with two ways to decompose the quartic vertex Tr BB^BB^ = tr B'^ BB^B with the help 
of auxiliary fields H and h, coupled respectively to BB'^ and B'^B and thus having the different sizes: 
N X N and n x n. 

As a word of precaution we remind only that for any finite n the Miwa transform p.lSP defines 
only an n-dimensional subset in the infinite-dimensional space of i-variables: when expressed through 
L the higher with k > n are actually algebraic functions of the lowest ti, . . . ,tn- Thus Z^it) in this 
context should be interpreted as a projective limit at n ^ oo. 



1.6 Genus expansion and the first multiresolvents 

Multiresolvents for Hermitian model are described in detail in the reference-paper |26j . Here we remind 
only the simplest of the relevant formulas. 

Multiresolvents are defined by eq. ijlOp and the first step is to rewrite Virasoro constraints as 
recurrent relations for p^P^'^\ Such recursive reformulation is possible only if the genus expansion of 
the free energy T = logZ is performed, T = X^^o S'^^"^-^*'^^ • explained in some detail in [26], this 
requirement picks up some rather special solutions to Virasoro constraints, and only such solutions 
possess well-defined multiresolvents and are associated with the bare spectral curves of finite generaH 
The bare spectral curve S is defined from non-linear equation for p^^\^^ - the starting point of the 
recursion. The next step provides p^^^'^\ which appears to be easily connected with Bergman kernel 
bi-differential on E [IS] . At each step of recursion there exists a certain arbitrariness in the choice of 
solutions, which is, however, absent for the Virasoro constraint p.5j) - crucial for this unambiguity is 
the form of the l.h.s. of ()I.5P : the fact that it is obtained by the shift of the time variable t2 — > t2 ~ 1/2. 
Shift is parameter which is allowed to stand in denominators when we build up a formal-series solution 
to Virasoro constraints. The shift of t2 is obviously associated with the integral (jl.ip and defines what 
is naturally called the Gaussian phase of Hermitian model. If other time or many times are shifted, 
then arbitrariness is unavoidable, see [49j for its full description. If partly fixed and parameterized by 
several arbitrary variables, it provides the family of Dijkgraaf-Vafa non-Gaussian partition functions 
[5Ol[5ll[l8]0 



Gaussian phase. The bare spectral curve for Gaussian phase of Hermitian model is 

^hg: y' = z^-^S (1.40) 
where S = gN and a few first Gaussian multiresolvents are: 

(o|2)/ ^ 1 f ziZ2 - 45 \ 

2{zi - Z2Y \y{zi)y[z2) J 

pf'W-^ (1.43) 

(0|3), X 2S{ziZ2 + Z2Z-i+ Z^Zi+AS) 

Ph {zi,z2,Z3) = 3/ ^ 3. X (1-44) 

y^[zi)y^{z2)y^{z3) 



^It deserves emphasizing that genus p in "genus expansion" refers to the genus of the fat-graph Feynman diagrams 
contributing to J^^K It has nothing to do with the genus of the bare spectral curve, throughout this text this genus 
wiU be only zero, while the genus of the full spectral curve (which defines the point of the Universal Grassmannian [JT] 
underlying the matrix-model r-function) is infinite. Relation between bare and full spectral curves is rather tricky and 
is not yet fully clarified in the literature. 

^Relation between generic genus-expansion-possessing solutions of and the Dijkgraaf-Vafa family is very similar 
to that between the "general" and "total" solutions to the Hamilton- Jacobi equation, see [521 sect. 7]. 
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s 



-^^—-^^——(ziZ2{bz\ + izfz2 + Szfz^ + 4zizf + 5z|)+ 
y\zi)yT[z2) V 

+45 {zf - 13ziZ2izl + ziZ2 + z|) + z^) + 16S'^{-zl + l3ziZ2 - zl) + 320S'- 



(2|i) 2lS{z^ + S) 



They are deduced from the recurrent relations 



y{zi)P^H'^\zi,Z2) = d. 



Zl - Z2 



y{zi)PH^^\zi) = Ph^^\zi,zi) 



.(0|2), 



y{zi)p%^^\zi,Z2,zz) = 2p'^^^\zi, Z2)p^H^\zi, Z^) + 
Smt^ . ^ J0|2)^^ ^ ^ J0|2). ^ (0|2), 

^ Ph \zi-,z2)-r 



„(0|2), 



J0|2), 



+5. 



pj^'^^ (21,2:3) - pff^^ (22,2:3) , a pij''^^(^i>^2) - pff^^ (22,23) 



- Z2 



Zl - Z'i 



y(2i)pi^'')(2i,22)=2pS"^(2i,22)piJI'^(2i) + pgl')(2i,2i,22) + a,,: 

y{z^)p^^\z^) = (f}^\z^)f ^p'f\z^,z^) 



Ph {^V - Ph K^V 



Zl - Z2 



(1.45) 
(1.46) 

(1.47) 
(1.48) 

(1.49) 

(1.50) 
(1.51) 



Non-Gaussian phases. For the sake of completeness we also give some formulas for non-Gaussian 
phases. If instead of the Gaussian shift i2 — *■ ^2 — ^ we apply ^ tk — Tk with W{z) = J2k=o "^kZ^ 
and rewrite the shifted Virasoro constraints in terms of the multiresolvents we get: 



where 



and 



W'{z)p{z) = p^(z) + f(z) + i\I{z)p(z) + P- lv'(z)p(z)] 



k=0 



dtk 



p{z) = il{z)g^ \ogZ 
f{z) = P+ [T^'(2)p(2)] = RH{z)gHogZ 

T^'(2)piJ™+^) (2, 21 , . . . , 2^) - /(^l-^^^) (2| 21 , . . . , 2^) 



{p-q\m2 + l) 



+E 

i=i '^^■^ 



— Pw ) • • • ) ^imi )PvF 

q mi+m2=m 

9 Pw \-^^ • • • i ■ ■ ■ i ^m) ~ Pw \zi,...,Zm) 



(2, 2j^, . . . , Zj ) + 



dzi 



2 - Zi 



+ V{z)p%-'^-^'\z,Z,,...,Zm) 



/^^^""^'Hzlzi, ...,Zm)= RHiz)p^P\"'\zi, . . . , 2™) 

Rh{z) = P+ \w'{z)V{z) 



d 



n—l n—a—1 

J2 E ia + b + 2)T,+b+2Z^^ 



(1.52) 

(1.53) 

(1.54) 
(1.55) 

(1.56) 

(1.57) 
(1.58) 



For further details we refer to 



and references therein. 
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2 Complex-matrix model Zc{t) 
2.1 Original integral representation 

Complex matrix model was originally defined as an integral over N x N complex matrices $ 

Tr ^ tk 



Zc{t) = I exp + V ^Tr {^^^f (1.59) 

JnxN \ 9 9 I 



where = (f)"^' Y{f,j=i d^^ij = n5=i dRe($,,)dIm($i,) = n5=i hd^^3d4r 
2.2 Eigenvalue representation 

One can express a complex matrix <I> through Hermitian H and unitary U matrices, 

= UH, $t = i^-f/t (1.60) 
and, further, through diagonal matrix D and two unitary matrices U and V: 

$ = UDV\ = yZ)[/^ (1.61) 

The norm of $ decomposes as 

Tr 5^5<^^ = TV {5Hf - Tr H^{U^6Uf + Tr [F, (5/i']([/U[/) = 
= Tr {6Df - Tr (C/tJC/)^!)^ _ Tr D^{V^6Vf + 2Tr {W6U)D{V^6V)D, 
so that the measure 



(1.62) 



N N 

= [dU] [dV] HiHf - H]f W HidHi (1.63) 

where Hi are the eigenvalues of matrix H (i.e. the entries of D). In particular, for = 1 we have 
$ = e^^i/ and d^^ = HdHdO. 

Comparing with p.2p . p.3p . one can see that for complex matrices the measure is actually the 
same as for Hermitian matrix , such that = UH'^U^. Since action in the model ()I.59P also 
respects this substitution, we obtain: 



~ J [dU] J exp (^-^Tr + ^^''^ ^(^^) 



(1.64) 



This integral looks just the same as (jl.ip for Hermitian matrix model, however, there is a difference 
in the integration contour: is not an arbitrary Hermitian matrix. Relation between p.ip and ()I.64p 
is like between an integral over entire real axis and over its positive ray: the answers are different and 
even invariance properties - and thus the Picard-Fucks equations (Ward identities) are not exactly 
the same. 



2.3 Virasoro constraints 

Since the Eq.(jL64]) means that the Virasoro constraints are the same as the "discrete Virasoro con- 
straints" for Hermitian matrix model. However, there are two differences. 
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First, L_i-constraint, associated with the shift 6{H'^) = e, is excluded, because it would correspond 
to a singular transform 6H ~ H^^. This exclusion can also be considered as a result of the above- 
mentioned change of integration contour: from entire real line in the case of Hermitian model to a 
positive ray < < oo in the case of ()I.64p . 

Second, the shift of time variables is tj^ = if^ — ^5k^2 in the Gaussian phase of Hermitian model, 
but it is rather tk = ik — 6k,i in the Gaussian phase of ()I.64p . 

The two changes together make the seemingly diminished set of Virasoro constraints 

LmZc = 0, m > 0, (1.65) 

with 
and 

-wr- = -2^c (1-67) 

fully exhaustive: despite the lack of these equations are enough for unambiguous recursive re- 
construction of all terms in the formal series Zc{t) for the Gaussian branch of the complex matrix 
model. 

2.4 Determinant representation and integrability 

The integrable properties of the complex matrix model are practically identical to those of the Her- 
mitian model. In particular, the partition function 

Zc ^ = deid'+^-^C (1.68) 
where the moment matrix is given a bit different integral as compared with the Hermitian case, 

Cij= [ dxexp I - VtfeX^ ) (1.69) 

This is still a Toda chain r-function, however, it corresponds to another solution to the hierarchy given 
by the Virasoro constraints ()I.65p . 

2.5 Kontsevich-Penner representation 

Like in the case of Hermitian model, the set of constraints p.65p has another matrix-integral solution 
[3T] . different from (|L59]) : 

Zc{t) = 7r^'-"'e-*^ ^''^ j exp (-tr (t)cp^ + tr r/V + tr (^"^r/ + Nil \og{^cp'^)) (fcf) (1.70) 

This time integral is over complex matrices (j), but their size n is, like in p.l7p . independent of N, 
which appears only as a parameter in the Penner term. For the sake of simplicity, we put here g = 1, 
the (/-dependence being easily restorable. The time variables are related to the external matrices rj 
and 7]'^ as 

tk = -7tr (wt)-^ k>l, 

^ (1.71) 

to = log(r?ry'^) 
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2.5.1 Proof I: Ward identities 

The Ward identity associated with the shifts 



d ^fd 

+ N 



dr]t, 



drj] 



+ T] 



tr, 



^ + e of the integration variable, is now 
j exp (-tr </>(/)^ + tr + tr (p^r] + iVtr log(#"^)) (fcp = (1.72) 



Therefore, one gets 
/ d 



d4r 



^_ 1 d__l ^ d 
dridrj^ 2^ drj 2^ drj^ 



Zc = 
(1.73) 



and, substituting Zc as a function of the Miwa variables ()I.7ip one reproduces (jL65]), see [3T| . 
2.5.2 Proof II: Orthogonal polynomials 

Let us put L = T]!]^ . Then, one can immediately repeat the calculation of s ]1.5.2l in order to obtain that 
(I1.59P is equal to p.24p and p.25j) . where the polynomials Pk{H) are now orthogonal with the weight 
exp(— x) on the positive real semi-axis. Such orthogonal polynomials are nothing but the Laguerre 
polynomials ^53], which have the following integral representation 



Pk{x) = {-ife- / dyy'e-yJo{2^) 
Jo 



(1.74) 



where Jq{x) is the zero order Bessel function. 

Now let us rewrite p.70p in the determinant form. This time in order to integrate over angular 
variables, we need to use instead of the Itzykson-Zuber formula the following very nice formula of 
integration over two unitary n x n matrices, [54] 



dU [ dFexp ( -tr 

Uin) Ju{n) V2 



A(a;^)A(y^) 



(1.75) 



where xf and y'j are the eigenvalues of A'^ A and BB'^ respectively, A and B being arbitrary n x n 
complex matrices. Using this formula and p.GSP and denoting eigenvalues of (j)(f)^ and r/r/^ through yi 
and Xi respectively, one obtains 



e ""^ / exp (-tr 



+ tr rycf) + tr (frrj + A^tr log 



~ e 



V ^ A(x) 



exp - ^ yi + iV ^ log yi Jo(2i/ 



det ^T'ix 



Xiyi) 



J, 



A(x) 



where 



$f)(x) =e-^' /"dyy*-iexp(-y + AflogyJo(2^)) 
Jo 

Comparing this with (jL74]) we ultimately identify (|1.59P and p.70p . 



(1.76) 



(1.77) 



2.5.3 Proof III: Faddeev-Popov trick 



Direct equivalence of the two integrals p.59p and p.70p can be proved by a somewhat tricky general- 
ization of Faddeev-Popov argument from [21j, which we applied in s ll.5l above. 
As before, we make the other choice of the sign in the Miwa transform. 



k > 1 



(1.78) 
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in order to deal with bosonic auxiliary fields. After the Miwa transform, the integral ()I.59P becomes 

r / °° 1 \ r p-Tr <I"I>t ,2^ 

(p9ll = / exp -Tr + y tr (r/r/t)-fcTr ($$t)fc ^ det(-W)A^ / ^ 

^ ^ -^ V j J det(^>^>t ®I -I® r/r/t) 

= det(-W)'^ / /e-Tr-f^t+tr Bt$$tB_Tr BWBt^2^^2^ ^ ^g^(_^^t)iv f g"^' '"''''''^'f, ^ 
J J J detArxAr(/ - -D-d"'') 

(1.79) 

Note that the last determinant is raised to the power N , this is because we integrate over A^^ complex- 
valued variables the relevant piece of the action is Y^fj^k=i ^ij^ik {^jk — J2a=i BjaBka)- 
At the same time integral p.70p is: 

^ det^x„((0 + r/)((/.t+r,t)) J J 7 det„xn(l + 6^6)" 

(1.80) 

Like B, b is rectangular N x n matrix, but the integrals are not literally equal as it was in the case of 
Hermitian model: one still needs to relate B and b. 
Let us begin with a few examples. 

Examples. N = n = 1 : In this case we can introduce new variables: p = B^B = |Sp and 
a = b'^b = Denoting also K = r/r/^^ = we obtain our two integrals in the form: 



N = n = l: (IL59]) = / while (OO]) = / ^ '^'^ (1.81) 

J 1 - p J 1 + a 



and integrands coincide because for p = we have dp = — ^^^^'^^ , while 1 — p = . 

n = 1, N arbitrary: In this case is a complex A-vector {Bi, . . . , -Bat), the A x A matrix BB^ 
has rank 1 and det(/ — BB^) = 1 — |-Bip — ... — |-BnP = ^~Pi~---~Pn = ^~ P+j so that 



and similarly 



e ^+<^^ da-i . . . dapj f ^+''^daA 



^ / e aai...aaN ^ / o+ e »a+ ^ (1.83) 

J 1 + (7+ J 1 + (7+ 

where we used the fact that the volume of a simplex pi > 0, . . . , pjy > 0, pi + . . . + p^ = p^ is 
proportional to p^~^ and similarly for cr's. Making the same transformation as above, = j^^, we 
see again that the integrals coincide. 

A = 1, n arbitrary: This time B and b are complex n-vectors. If we perform an SU{n) rotation to 
diagonalize r]r]^ — > diag(A'i, . . . ,Kn), then both integrals p.SOP and (|I.70P still contain B and b only 
in the form of squared modules pa = |^aP and = \ba\^: 

J 1-/9+ J (1 + (T+)" 

The integrals are related by our usual change of variables pa = , only the measure transform gets 
a little trickier: 



f\''a=ldpa = /\a=l 



dag aada+ \ _ A2=id(Ta Ea=i o'a ^^"=1 ddg _ 1 

l + a+ {l + a+yj (l + a+)" (l + c7+)"+i (1 + a+)-+i ^^=1 

(1.85) 

Substituting also 1 — /9+ = (1 — cr+) we see that the two integrals are in fact the same. 
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2.6 Genus expansions and the first multiresolvents 
Gaussian phase. The first resolvent [26] in this case is 



and the bare spectral curve 

Sc : y^ = z{z - 45) (1.87) 

A few next multiresolvents are: 

(0|2). ^ 1 ( ZiZ2-2S{zi+ Z2) \ 

py(.„z„.3) = 2 3^^yf; (1.90) 

yc(2;i)''yc(^;2)'^yc(^;3)'' 

v2^2 c2 



, , ^ ~2^2q2 

(i|2)^^. _ ziz^b ^ ^ 



/9^''"^(Z1,Z2) = (3Z2^? + ^Zlzi + 2z|z2 



yiy2 ^ (1.91) 

-2{zl + nzfz2 + nzizl + zl)S + 8(3z? + l9ziZ2 + 2,zl)S'^ - 96(^1 + Z2)S^ + 1285^ 

^(2|i)^^^ _ (95^ - SzS + 8z^)g^z3 
2/(2;) 

They are deduced from the recurrent relations 



pgi^)(,) = ^ (1.92) 



^pgl^)(zi,Z2) = (1.93) 

Zl Zi Zl- Z2 

y-^p^\'\z)=p^^'\z,z) (IM) 
^pgl^)(z„Z2,.3) = 2p^^^'\z„Z2)p?'\z^,z,)+ 

1 (0|2)/ X (0|2). X , (0|2). X (0|2), X (^-^^^ 

, 1 _f) ^iPc (^l;^3)-^2Pc (^2,^:3) , 1 o ^iPc (^1^^2) - ^^3^(7 (^2,^:3) 

Zl Zl - Z2 Zl Zl - Z3 

(1|2)/ X o (0|2)^ ^ ^ , (0|3). \ , 1 ^IPC ^^(^1) - ^2Pc ^^(^2) 

/'C {ZI^Z2) =2p}j' '{zi,Z2)p'c' '[Zi)+p]^' '{zi,Zi,Z2) + —0^2 

Zl Zl Zl - Z2 

(1.96) 

^pg'^^(^i) = (pg'^H^i))%pg'^^(^i,^i) (1-97) 

Non- Gaussian phases. A generic phase of the complex matrix model is given by first several time 
variables shifted, ~^ + tk- Then, for a generic polynomial potential W{z) = J2k=o'^kz'' the 
Virasoro constraints for the complex-matrix model look like 

W'iz)p{z) = p\z) + f{z) + 9^V(z)p{z) + [zv'{z)p{z)] (1.98) 

where 

f{z) = ip+ [zW'{z)p{z)] (1.99) 
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Higher resolvents can be extracted from the equations 

— ^ 2:^^ , . . . , Zj^^ ^' ^(z, Zj^, . . . ,2;^^^)+ 

5 m,i+m2=m 



where 

ic 

with 



. . . , z^) = Rc{,z)p^\^\z,, ...,zm) 



n—l n—a—1 



Rc{z) = -Pt zW'{z)V{z) = - E E {a + h + 2)Ta+,+2Z^^ 

^ a=-l 6=0 



Quadratic equation for simplest resolvent leads to the answer 
where 

yc{zf = z{w'^-mc{z)F§'^) 
For the Gaussian complex model = 5fc,i, so W{z) = z, Rc{z) = ~jgf^ 
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3 Kontsevich model Zk{t) 



3.1 Original integral representation 



Kontsevich model was originally defined in |22] as a generating function of topological indices of the 
moduli space of Riemann surfaces. M. Kontsevich represented this generating function in form of 
the now-famous matrix integral over auxiliary nxn dimensional Hermitian matrices (one can easily 
introduce into this integral the parameter g similarly to p.l7p . [1] however, for the sake of simplicity, 
we put here g = 1): 

where time- variables are Miwa-transformed: 

Tfc = \ti L-^ (1.106) 
k 

If expressed through the r- variables, Zk{t) is actually independent of auxiliary parameter n. This 
model can be further generalized to Generalized Kontsevich Model (GKM) [23j , 

J g-l/(L,h)^^ 

Zgkm = -p^ — ^^ (1.107) 

where 

U{L, h) = tr [V{L + h)- V{L) - V{L)h] (1.108) 

and 

U2{L, h) = lim ^U{L, eh) (1.109) 

is an /i^-term in U . W is here an arbitrary (power series) potential and Zgkm is a function of the 
same (VF-independent) Miwa transform (jLToHj). Many properties of GKM are in fact independent of 
the choice of V{h). 

In fact, one can also consider the matrix model 0.107^ with a different normalization as a function 
of time variables 

% = itr (1. 110) 

where the matrix A = V'{L) enters in positive powers. This function is called the character phase 
and is considered in detail in [20]. In this paper we restrict ourselves with the Kontsevich phase only, 
where the time variables are given by p.l06p . 

Note also that it is often convenient to fix V(/i) to be a polynomial of h (polynomial Kontsevich 
model, [23j) or that of (antipolynomial Kontsevich model [20j). In this section we consider only 
the polynomial case, leaving the antipolynomial one until the next section (where it emerges within 
the context of the unitary matrix model). 

3.2 Eigenvalue representation 

Shifting the integration variable h ^ h — L one obtains that 

Zgkm ^ I dh exp {-tr V{h) - tr V'{L)h) (I.lll) 

J nxn 

Now using the Itzykson-Zuber formula and ()I.3p . one can perform integration over the angular variables 
in this integral: 

/ dh exp (-tr V{h) - tr V {L)h) ~ ^^^^i^'^^i) (1.112) 

Jnxn ^(-^j 

where Aj are the eigenvalues of the matrix V'{L) and 

Fi{X) = J dxx'-^ exp (-V(x) + Ax) (1.113) 
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3.3 Virasoro constraints 



Straightforward Ward identities for Zk are as previously associated with the shift h 
integration variable h: 



/i + e of 



d 



dLl 



dh exp 



--tr h!^ +tr L^h 
3 







(1.114) 



Now one should take into account the normalization factor and come to the r-variables. Conversion 
to the r-variables is highly non-trivial, it was first performed in [18] and leads to the celebrated result 

[nmi]: 



+On+l,0 • —r + "n,0 ■ — 

4 lb 



a + b— 2n " " 

a,6>0 ; a,b odd 

d 



+ 



(1.115) 



This proved equivalence of Witten's topological 2d gravity [12j to Zk and - since Zk is trivially a KP 
tau-function - proved that partition function of 2d gravity is indeed a tau-function (as anticipated in 
[55]). Analogous conversion to T- variables of the Ward identities for Zgkm is even more sophisticated 
and give rise to Ty-constraints (or M^-constraints in the character phase) [191 120j . 



3.4 Determinant representation and integrability 

Now one can take into account all the normalization factors and further transform this determinant 
(after quite tedious calculation, [23]) to 

ZOKM = ^^^^ (1.116) 

where 

UL) = yV%^e^(^)-^^'(^)Fi(V'(L)) (1.117) 

Formula p.ll6P if true for any number of Miwa variables (size of the determinant) fixes a KP hierarchy 
r-function [THl [24] that depends on times 

^k = \T.L^^ (1.118) 

i 

provided the asymptotics of (piiL) at large L is 

0,(L)''~°°L*-i(l + O(l/L)) (1.119) 

In particular, (1I.119|1 guarantees that Zgkm is a function of variables (1I.118|) and does not depend 
on their number. 

Now if one takes the monomial potential V{h) = h^^^ (the case of a polynomial potential of degree 
p + 1 describes a hierarchy equivalent to the monomial case, see details in [56]), the partition function 
Zgkm is a r-function of the p-reduced KdV hierarchy, which does not depend on times r^^ for all k. 
In particular, the Kontsevich partition function p.l05p is a KdV r-function depending only on odd 
times T2k+i- A concrete solution of the KdV hierarchy is fixed by the Virasoro constraints ()I.115P (in 
fact, it is enough to use only the lowest constraint in addition to the KdV hierarchy equations in order 
to fix the partition function unambiguously). 

Note that one could starts from the Virasoro constraints 0.115^ instead of the matrix integral. 
Then, there are much more solutions, the KdV one corresponding only to distinguished solutions of 
the Dijkgraaf-Vafa type [26', 2nd paper]. 
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Note that one can easily continue the (Generalized) Kontsevich matrix integral to the whole Toda 
lattice hierarchy adding to U{L, h) (but not to C/2(L, h)) the term 

AJ7(L,/i) = AV(L + /i)-AV(L), AV(/i) = Klog/i - ^fjk/i-*= (1120) 

k 

Here H is the zeroth (discrete) Toda time and fk are the negative Toda times. In the special case of 
quadratic potential V{h) = h"^ this matrix integral reduces to the Toda chain, as we observed in s.2. 

3.5 Kontsevich- Penner representation 

Of course, Kontsevich model is already in the Kontsevich form. No parameter is obligatory present 
and no Penner term is needed (until one wants to deal with the (Generalized) Kontsevich integral as 
with the Toda lattice hierarchy). 



3.6 Genus expansion and the first multiresolvents 

Generic phase. Similarly to the Hermitian and complex matrix models, a generic phase of the 
Kontsevich model is given by first several time variables shifted T2k+i — T2ifc+i + Tfe. In this case. 



v(^) = E 

k= 



^ Z^+^ dT2k+l 



n+l 



W\z) = J2[k + i,)TkZ 



k=0 



1 



oo . 
k=0 ^ 

v'{z)p{z) 



W'{z)p{z) = p\z) + fK{z) + g'V{z)p{z) + zP, 



, 9' , (ro-To)^ 
16z 16 



fK{z) = ZP+ 



W'{z)p{z) 



Rk{z) log Z 



n+l fc— 2 ^ -. ^ 

RK{z) = -J2J2[k + l^]Tkz'-^ 

k=2 m=0 



dTrr 



i=l 



p{z)=g^V{z) log Z 

W'(z)pi^™+^) {Z,Z,,..., Zm) - flP'^^'Hzlz^, ...,Z^) = 

E E Pw \^i^hi ■ ■ ■ i^im'i)Pw ''^^ '*(^) ^ji ) ■ ■ ■ ) ^jma ) + 

q mi+m2=m 



1_\ ZiP^^ '{z,Zi,...,Zi,... , Z 

^KdZi 2Zi 



{p\m) / N 
Zpw l^li---i^mj , ^/ X (p-l|m+l) 



(1.121) 

(1.122) 

(1.123) 

(1.124) 
(1.125) 

(1.126) 
(1.127) 



Z - Zi 



+ V{z)p)^-'^"'^"{z,Z,,...,Zm) + 



16z 



. ^ ( SmflTf 6ni,lTi 6m,2 

+ 8^ + 8^2 



(1.128) 
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Gaussian case. In this case, W'(z) = z + '^, /('^l™) = 

yoii)(,)_ ^ + f -^^^+^ 

(012)/ X 1 f + Z2 + 2Tq)ziZ2 . \ 

V a2;2 2^2 y z\- Z2 8^2 



p(0|^)(zi,Z2,Z3) = 



z\z2z\ 



y(zi)3y(z2)3y(z3)3 



y(zi)p(0|3)(2i,Z2,Z3) 



dz2 1Z2 J Zi- Z2 

8ziZ2 



(_d ^\ zsp(''\^\z,,Z2)-zip(^\^\z2,zs) ^ J_ 

\dz3 2z-i) ZI-Z2, Szi; 



y(.)p(iW(z)=p(0|2)(z,.) + ^ 



,(i|i)(^) 



y(zi)p(i|2)(zi,Z2) = 2p(°|2)(zi,Z2)p(^l^)(zi) +p(°l^)(zi,zi,Z2) + 



lQy{zf 

d 1 \ Z2P^^\^\zi) - Zip 



^ 52;2 22;2 y zi- Z2 

y(zi)p(2|i)(z) = (pWi)(^))Vp(i|2)(^,^) 

(211), X ^^|5(zi+To)2 + 3(zi+To)(z2+ro) + 5(z2+To)^ 
^(-i'-2) = ^ y{ziyy{z2y 



254y(z)ii 
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4 BGW model Zxir) = Zbgw{t) 

4.1 Original integral representation 

Brezin-Gross-Witten (BGW) model is defined as a generating function for all correlators of unitary 
matrices with Haar measure \dU\. 

Zbgw = I [dU\ exp fXr J^U + Tr JU^) (1.140) 

JNxN ^ ' 

The integral actually depends only on eigenvalues of Hermitian matrix M = Jj\ i.e. on the time- 
variables of the form = Tr (JJ^)^. 

Haar measure [dU] for unitary matrices is non-linear, it can be reduced to a flat measures in different 
ways. One possibility is to express U through Hermitian matrices, U = ^5'7j , which defines [dU] 

as the flat Hermitian measure dH = Ylfj=i dHij with additional Jacobian factor, [dU] = J'{H)dH, 
J = det{l + H^)-^. Another possibility p0l [2T] is to impose the constraints on the complex matrices: 

[dU] = f d2$5($$t -I)= [ dhe-^^ ^ I d^^e^ (1.141) 
J Jnxn Jnxn 

For certain actions the integral over d^ can be explicitly taken and this gives rise to reformulation of 
original unitary-matrix model. 

4.2 Eigenvalue representation 

Since technically the most simple way to obtain eigenvalue representations is to start with the 
Kontsevich-Penner representation of the BGW model, we first consider this representation. 

4.3 Kontsevich-Penner representation 

In variance with all other Kontsevich-Penner representations, that of the BGW model connects the 
two integrals over the two matrices (unitary and Hermitian ones) of the same size, |20] : 

Zbgw = I dhexp( - tr h'^ + ti Mh- Ntr log h) (1.142) 
Jnxn ^ ' 

This makes theory of the BGW model somewhat harder and one sometimes embeds it into the universal 
BGW model [20] with an arbitrary coefficient in front of the logarithmic term. However, in order to 
make contact with the BGW model (|I.140p one ultimately has to put this coefficient equal to —N. 

4.3.1 Proof 1: Ward identities 

The simplest Ward identity for Zbgw has the form 

Zbgw{J, J^) = I- Zbgw{J, J^)- (1.143) 



djl dJtr 



or, equivalently, [2D] 

At the same time, integral (II.142P satisfies the equation 



dMtr dMtr dMtr \dMtrJ \dMtr 

AtM = N and V'{h) = l/h p.l44p and p.l45p coincide, which establishes ()I.142p . 



^^gtr (Mh-Af log h+Vih)) ^ Q 

NxN 

(1.145) 
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4.3.2 Proof II: Faddeev-Popov trick 

Another simple way to derive the Kontsevich-Penner representation for the BGW model is to use the 
trick (|LT4T1) . [H]. Indeed, 

Zbgw = J dhe-^^' ^ J d^<^ exp (tr + Tr + tr J^^^ = 

= / d/i exp f - tr /i + tr M /h - Nti log h) ^^=^ (I 146) 

Jnxn ^ ' 

The BGW model has two phases |20j: the Kontsevich phase, where partition function is expanded 
in negative powers of M and the character phase where expansion goes in positive powers of M. Below 
we describe them separately. 

4.4 On direct relation between Zq and Zbgw 

In the Kontsevich-Penner form p.70P the complex matrix model looks somewhat similar to original 
form (II.140[) of the BGW model. From (II.70|) one obtains (representing cf) = HU'^ and 0^ = UH): 

This tricky formula is the best direct relation known at present. A more transparent relation is still 
lacking. 

4.5 Character phase 

In this phase the BGW partition function is considered as a function of the variables 

Tfc = ^tr M'' (1.148) 
k 

and one has to consider the Universal BGW model, i.e. the Kontsevich integral ()I.142p with an 
arbitrary coefficient of the logarithm, which is a free parameter and not the size of the unitary matrix. 

4.5.1 Virasoro constraints 

Performing the change of variables in (II.145jl from M to 7^, one can directly obtain the Virasoro 
constraints satisfied by the BGW partition function Z^^y^^ in the character phase: 

Lni{N,T)Z^Qy^, = 5m,iZ^(j^, m>l 

d °° d (1.149) 

Therefore, the Ward identity (and its solutions) depends on the size of matrix N. This means that 
the integral p.l42p is not just a function of variables T, but also depends on N. The way out is to 
consider the Universal BGW model given by the integral 

ZuBGW = dhexp( - tr /i"^ + tr Mh- Afti log h) (1.150) 

JnxN ^ ' 

Then, the Virasoro constraints become 

L„(2Af - AA, T)Z^BGW = ^m^Z^BGW^ ^ > 1 (1-151) 

and choosing J\f = 2N — H, one arrives at the partition function Zbgw that does not depend on 
(but only on the parameter N) though the integrand in (|I.150P does! 
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4.5.2 Determinant representation and integrability 

One can easily integrate over tlie angular variables in ()I.150p as in the previous section to obtain 

detFi(M,) 



where 



ZUBGW = (1-152) 

Fi{M) = j dhh'-^ exp ( - + MH - M\og /i) = 2m (2\/m)'^~* lAT-i (2\/m) (1.153) 

where Ik{z) are the modified Bessel functions. After some work [20], this formula can be recast to the 
form p.ll6P with the asymptotics (jl.llOp . where L = 1/M and 

T{2N -Af-2i + 2) /^\^^~^'\ f2\ 

— : — : — ^ J^J_^/'/9_,• — = 



Again in order to make these functions independent of A^, one has to choose M = 2N — b^. At the 
same time, this proves that, under such a choice, Z^^q^^ is a r-function of the KP hierarchy. 

4.6 Kontsevich phase 

In the Kontsevich phase the unitary matrix integral is considered as a function of variables 



Tfc = -itr M-'' (1.155) 
k 

This time the integral p.l42p does not depend on A^ provided it is properly normalized: 

-2 tr 



^BGW ~ ^ 



det (Ml/2 M + M MV2) 

\ i^:^ ^^^^ ^'-'"'^ 



One can check by a direct (quite involved) calculation [20] that Z^^^^^ depends only on odd times 

4.6.1 Virasoro constraints 

Using the Ward identity p.l44p one can now make the change to variables (|I.155p to obtain the 
Virasoro constraints satisfied by Z^Q^y [331 ED] : 

^mZsGW = 0' rn>0 

^o^fc ^^^^k-2m 2 ^J^^ dTadn dT2,n+l 16 

a + fe— 2m 

4.6.2 Determinant representation and integrability 

In the Kontsevich phase, performing integration over the angular variables in p.l42p and taking into 
account the normalization factor, one obtains [20] the determinant representation p.ll6P with the 
asymptotics p.llOp . where 



(M) = 2,/^e~'^^^ M'-^-^/^I,^N~ii2M) (1.158) 



This means that the partition function Z'^Qy^ of the unitary matrix model is a r-function of the KP 
hierarchy. Moreover, as it was already noted, it does not depend on odd times and is, in fact, a 
r-function of the KdV hierarchy. 
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4.7 Genus expansion and the first mult iresol vents 

Generic phase. As before, we shift first time variables T2k+i T2k+i + 

n+l 



^ / 1 \ 



and obtain the loop equations 



W\z)p{z) = p\z) + fBGw{z) + g^V{z)p{z) + P" + ^ 

/i^(^) = P+ [W'{z)p{z)] = g^RK{z)\ogZ 

n+l fc-1 /IX o 



,2 



27 9T„ 



=52V(z)l0gZ 



16z 



1.159) 

1.160) 
1.161) 

1.162) 
1.163) 
1.164) 
1.165) 



W'{z)p^tr^'\z, .1, . . . , .^) - M'^^'\zW. ...,z^) = 

— ^(^'^il' ■ ■ ■ )^2mi)Pl^ ^' '*(^5^jn • • • J^jma)^ 

g mi+m2=m 



(1.166) 



Gaussian phase. In this case, W\z) = 1, /C^l™) = 

/9(°|i)(z) = (1.167) 

/5(°l')(^i,^2) = (1.168) 

p(°l'H^i,^2,^3) = (1.169) 

p('l')(^) = T^ (1-170) 
Ibz 

p«l>(2) = (p<l|l>(2))' + pW2)(2,2) (1.172) 

= ^ (U74) 
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Free energy in the Gaussian case. As a direct corollary of Virasoro constraints (|I.157p . one can 
calculate the free energy expansion in the parameter g, log Zbgw = Yl,k=o9'^^~'^-^BGw'- 

^BGW = -'^/8 ln(ro-2) 

^(2) _ 9 Ti 



32 (ro - 2)3 (1.175) 



(3) ^ _225 T2 567 Ti^ 

■^BGW g4 _ 2)5 + 64 ^T-Q - 2)6 

(4) _ 55125 rg 388125 T2n _ 64989 

•^BGiy - -^YT (ro - 2)7 + ^12" (tq - 2)8 " ~64~ (tq - 2)9 

In general -^jjc^i/ is a polynomial 

^BGW = X! Cki,...,k„,Qkl ■ ■ -Qkm (1.176) 

fcl+... + fcm=p— 1 

of the variables = ^^^_^^2fc+i • This is the best illustration of drastic simplicity of the BGW 

partition function as compared to the Kontsevich and Hermitian cases, where all J-^^^ are sophisticated 
transcendental functions, and are simplified only in terms of moment variables. One may say in the 
BGW case the moment variables are extremely simple. 
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Part II 

Decomposition formulas 

1 The idea of decomposition formulas |[1| 

The key observation is that multiresolvents - if defined according to the rule (|lUp - are polydifi'erentials 
on the bare spectral curve S, intimately related to the U{1) current J^{z) on S, with prescribed 
singularities: usually they are allowed at some fixed points {punctures) on S. In this approach the 
Virasoro constraints on partition function are written as 

P_ [j\z)) Z = i JC{z, z') (j^{z)) Z = (III) 

with a certain kernel IC{z,z'), made out of the free-field Green function on E. The current is also 
"shifted": J^{z) — > J{z) + /S.J{z) and partition function Z depends on the choice of: 

• the complex curve (Riemann surface) S, 

• the Green function IC{z,z'), i.e. projection operator "P., 

• the punctures on S and associated loop operator J'{z), 

• the local coordinates in the vicinity of the punctures, 

• the involution of the curve with punctures and loop operator, 

• the shift AJ'{z) on S, 

• the contour C which separates two sets of punctures. 

If contour C goes around an isolated puncture, Z is actually defined by its infinitesimal vicinity 
and depends on behavior (the type of singularity) of J{z) at this particular puncture. Coordinate 
dependence is reduced to the action of a unitary operator (Bogoliubov transform, and exponential 
of bilinear function of J^) on Z. Types of singularities and associated Z's can be classified, and our 
quartet Zh, Zc, Zk and Zbgw are the lowest members of this classification. The former two are 
associated with a puncture at regular point of S, while the latter two - with that at a second-order 
ramification point. Zc and Zbgw differ from Zh and Zk by the choice of projection operator "P., 
i.e. the kernel JC{z,z'). 

If contour C is moved away from the vicinity of the puncture, it can be decomposed into contours 
encircling all other punctures: this provides relations between Z's of different types, associated with 
different punctures. If S has handles or boundaries, there will be additional contributions, associated 
with non-contractible contours - the corresponding elementary partition functions are not yet identified 
and investigated - this seems to be a very interesting problem of its own. 

In what follows we present the two simplest examples of this procedure, both associated with 
S, represented as a double-covering of the Riemann sphere with two ramification points. Such S 
is of course also a Riemann sphere, however, representing it as a double-covering provides a simple 
description of behavior, which we allow J^{z) to have at the two ramification points. The other pair 
of punctures are chosen at preimages of a regular points {z = oo± in what follows). After that, 
depending on the choice of projection operator V- we obtain either a relation between Zh and the 
two Kontsevich models, Zh = Ukk(^Zk ^ Zk^, or between Zc and the pair: Kontsevich model and 

BGW model, Zc = Ukbgw(Zk ® Zbgw^- These both examples were already described in pj, but 
here we provide a more targeted and, hopefully, more clear presentation of the subject. Some mistakes 
of original version are also corrected, in the case of discrepancies from [T] the present version should 
be trusted more. 

2 The basic currents, shifts and projection operators 

These are the data, defining the standard Virasoro constraints ^ and ^ and thus the four models, 
discussed in the section |T] above. All the four are defined in vicinity of a particular puncture and do 
not depend on the global properties of the bare spectral curve S. 
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2.1 Hermitian current 



This one is used in the definition of Q and thus of partition functions Zn and Zc- 

fc=0 



(11.2) 



With this current one can immediately associate a bi-differential fj{z,z') = J{z)J{z')— : j{z)J{z') : 
where the normal ordering means all placed to the left of all t-derivatives. It is related to the central 
extension U{1) and is equal to 

/ 1 2 \ 2 dzdiZ 



fHiz,z'\g')=g' 



(11.3) 



2{z-z')^ 

This bi-differential will play an important role in comparison of global and local currents and, therefore, 
in construction of conjugation operators in the next subsections [3] and HI 

The further difference between various partition functions comes from different choices of the shift 
functions W{z) [26\ and projector operators [Tj, 



p 

J. 77- 



E 



ttkdz 



E 

k=m 



akdz 



(11.4) 



The two correlated choices lead to the two simplest models, associated with (jTL2]): to Zh an Zq- 

1. Gaussian Hermitian model [26j 

This model corresponds to the shift 



A ^ / N zdz 
^Jh{z) = - — 

Partition function is completely fixed by Virasoro constraints 

fH{z)ZH = 



where 



Lr, 



n=— 1 

52 



E^i^fc 

k=l 



2 



d 



dt 



k+n 



+^^E 



fc=0 



dtkdtn-k 



with 



^Znit) = ^Znit) 
oto g^ 



(11.5) 



(11.6) 



(11.7) 



(II. 



Given (jll.Sp . the choice of P_i from all the Pm is distinguished: with this choice only partition 
function is unambiguously defined by ()II.6p . There are interesting situations, when the choice of 
Pm is not adjusted to the shift in this way: the best known example is provided by Dijkgraaf- 
Vafa partition functions \50\ I51j. where projector is the same P_i as in Gaussian model, but the 
shift AJh{z) = dW{z) is generated by polynomial VF(2:) of degree higher than two. 

2. Gaussian complex model |30] 

This model corresponds to the shift 



A ? / N dz 
AJc(^) = - Y 

Thus partition function is completely fixed by Virasoro constraints 

fc{z)Zc = 



(II.9) 



(11.10) 
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where projector is taken to be Pq - again, to guarantee the uniqueness of the solution to (jll.lOp . 
- and 

{dzf 



fc{z) = Po [. {Jh{z) + AJc{z)f :J = 5' E 7^2 

d 



k=l 



dtk+n 



n=0 
Q2 



k=0 



dtkdt, 



n—k 



with 



(11.11) 



(11.12) 



2.2 Kontsevich current 

This one is used in the definition of ([6]) and thus of partition functions Zk and Zbgw^ 



2 _d_ 



(11.13) 



is - up to traditional but unimportant change of time- variables - the even part of the current pi.2p . 
However, associated central term bi-differential looks more sophisticated (being the symmetric part 
of the bi-differential): 



jK[t,,t. \g ) - 9 



(11.14) 



4(^2 _ ^/2)2 

The simplest partition functions, associated with this current, are Kontsevich r-function and BGW 
model. 



1. Kontsevich r-function [22], [23] 

This time the shift is 



K 



(11.15) 



and the relevant projector is P-2'- 



= E ( ^ + 

A:=0 



TK{i) = P-2 
1 



n=-l ^ 



Tk 



d 



n-l 



k+n 



+^^E 



92 



fc=0 



5Tfc9r„„i_fc 16 16^2 



(11.16) 



Then Zx is uniquely defined by 

2. Brezin Gross Witten model |20| 

Now the shift is 



Tk{OZk = 



AJbgw 



cdj 



and 



n=0 ^ 



d 



n-l 



92 



+ 



16 



(11.17) 



(11.18) 



(11.19) 



with projector Pq unambiguously specify Zbgw by 

TsGwiO^BGW = 



(11.20) 
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3 Decomposition relation Zn — > ® Zk 

Now we can select a bare spectral curve S: 

y^ = z^-a' (11.21) 

select the punctures: at z = ±a and z = (X)±, and select the global current by allowing specific 
singularities at punctures: 

J{z\a^) = ^ ("fc + (^fc + zB^)y'^-^dz + g'iCk + zD^)^^ (11.22) 

A:=0 \ y 

Bi-differential for this current 

_ 2 (2:^' - a^)dzdz' 

- ' 2[z - z'Yy{z)y{z') ^"-''^ 

is defined by commutation relations 

2 d k+l d ^ d ^ ^ 

dAk A; + 1 dAk+i dBk ^ ' 

At punctures it is equivalent to the bi-differentials of the basic currents: 

fj{z,z'Y-^^ fn[z,z') 
fjiz,z') '^^^UfKiC±,e±) 



(11.25) 



where are some local coordinates in the vicinity of ramification points a and —a, defined respectively 
by 



00 

z = a ' 

k=l 

and 



+ E »te+ (11.26) 



00 

z = -a + J2 «fc ^-'^ (11.27) 

k=l 



The current pi.22p itself is equivalent to the currents from sl2J 

J{z) ^-^^^ Jh{z) 
J{z) ^-^'^ 2J^(0 

Time-variables in parametrization of the global current are related to local time as follows: 

k I z 

d 1 



(11.28) 



(11.29) 



z^J 



Tk 



dtk g"^ 

J 



2k + lJ ^2fc+i 



(11.30) 



Global current is related to local currents by conjugation operators. Conjugation operator at infinity 
is 

Uh = -A / {fj{z,z')-fH{z,z'))aH{z)(lH{z) = 



^ i i Ph^^\z,z')v{z)v{z') 
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where pn is a bi-differential counterpart of the two-point function of Gaussian Hermitian model 

= 1 {fj{z,z') - Mz,z')) = ^ ( TVA - ^] ("-32) 
r 2(zi - Z2y \y[zi)y{z2) J 

At ramification points the conjugation operator is as follows 

^H = ^J2 f f {fAz,z')-45,,fKiC^,C,))^K{mKiCJ) (11.33) 

To establish required Virasoro constraints one should shift the global current 

which leads to a shift of the conjugation operators: 

UH^Uu + ^j ^-^^hH{z)dz = U+^ f p'§^^\z)v{z)dz 

Vu^Vn + j, ^^^^ j + f^_^^ ^ j o^(e-) j 

(11.35) 

Then the projector 

with contour C encircles the segment ramification points ita on the spectral curve (but not the point 
z\ so that always \z\ > \z'\) do the job: since 

{z'f 



1 {^r 

(z - z')dz' f-'^ z^+^dz' 

^ ' k>0 

it picks up the terms with n > — 1 in infinity and since 



(11.37) 



1 _ (g^)2fc-l 

,A ~ ne2k+2^e' (11.38) 



{z-z')dz' ^^2e''+^d(' 



it picks up the terms with n > —1 for ramification points. 
After all we get the decomposition formul^l 

Znit) = e''"e''^ZK{T+)ZK{T-) (11.39) 

4 Decomposition relation Zc — > Zk ® Zbgw 

This decomposition formula was the topic of s.8 of ref.[I], however, it is described there in a too 
sketchy and partly misleading form. Thus we provide here a more detailed and careful presentationj^ 

^Actually, as it was already indicated in PP, we get a whole family of such formulas, with infinite set of free parameters 
given by coefficients aj: in (|II.26|) . (|II.27|I . 

® In [1] we considered decomposition formula for the complex model, starting from the same spectral curve pi.21|l as 
for the Gaussian Hermitian matrix model 

yl = z^- 45 (11.40) 
with additional puncture in z = 0±. However, the global current which we introduced was singular at 0±. Actually, in 
notations of [T], the proper global current should be defined on the curve 

yl^z^z^-'^S) (11.41) 

It is more natural to consider instead the current on 

yc = e(C-45) (11.42) 

of which the previous one is a double covering z = as we do in the present text. 
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The bare spectral curve is 

yl = z{z - 45) (11.43) 

and the four punctures are chosen at z = 0,45, oo±. Accordingly on this curve we define the global 
current 

J{z) = Y^(k + -){Ak + zBk) y^c^-'dz + g^^^ (Ck + zDk) (11.44) 

fc=o ^ ^ yc 

with commutation relations 



(11.45) 



dBk dAk 
and the global bi-differential 



fjiz,z') = g 



^{zz' -2S{z + z'))dzdz' 



2{z - z'fycy'c 
_ , 2yyS + {zz' - 2S{z + z') + ^S^){yl + yg) ^ , 
^ 2{yl - y'SYycy'c 

At punctures this bi-differential is equivalent to the following canonical bi-differentials from sl2l 



(11.46) 



fj{z,z')'^'UfK{z,z') (11.47) 

fjiz,z') '^UfKiz,z') 

The current has the following behavior: 



(11.48) 



Jiz] 


1 ~ Jh{z 


J{z) 


''^^ 2Ji^(^4 


J{z) 


2j^(e- 



with ^_ - some local coordinates in the vicinities of 45, respectively: 



= 45 + ^ alt 

k=l 

oo 
k=l 



2k 



(11.49) 



Time- variables of the local currents are expressed through those of the global one in the same way as 
in sl3l 



d ^ I k^ 
dtk g^ J 

2 r J 



(11.50) 



(11.51) 
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Global current is related to local currents through conjugation operators. Conjugation operator at 
infinity 

UC = -A <f {fjiz,z')-fHiz,z'))nHiz)nHiz') = 

9 Joe Joe ^jj ,.^) 

— i i p'§^'^\z,z')v{z)viz') 



2o2 

where pc is a bi-differential counterpart of the two-point function of Gaussian Hermitian model 

= 1 {fjiz,z') - fHiz,z')) (11.53) 

At ramification points the conjugation operator looks as follows 

2 

The shift of the current 



^C = ^Y1 / / {fj{z,z') -4SijfK{^i,^j))^K{^i)^K{^j) (11.54) 



(11.56) 



j^j-y^ (11.55) 

corresponds to the shift of conjugation operators 

Un^Uu + ^j ^^^nH{z)dz = U + \j p^^\^\z)v{z)dz 

The difference from the case of Hermitian model is that now we should get 

Ln, n > at oo, 

£„, n > at a, (11.57) 
J~-n, n > —1 at — a 

Thus this time the proper projector is 

^'j-n ■■ iJiz') + A J(z'))' : (11.58) 



Tc {z — z')dz' 

an we finally obtain the decomposition formula for complex model: 

Zcit) = e''^e''^ZK{T+)ZK{T-) = e'''^ e""^ ZK{T+)ZBGw{r-) (11.59) 

Conclusion 

In this paper we demonstrated that decomposition formula Zh — > Zk ® Zk of partition function 
for Gaussian Hermitian model into two cubic Kontsevich models has as its closest analogue another 
decomposition: Zq Zk ® Zbgw of the Gaussian complex model into the cubic Kontsevich and 
Brezin-Gross-Witten models. Thus all the four models are indeed the very close relatives, though this 
is not quite so obvious from their original matrix- integral representations. This paper is therefore 
an important outcome and summary of many different approaches, worked out during the years of 
development of matrix-model theory. It brings us one-step closer to providing a unified look at the 
whole variety of eigenvalue models and building up the M-theory of matrix models, suggested in [2]. 

Technically it adds to content of an identification of partition function, denoted there by Zk, 
with that of the very important BGW model - the generating function of all unitary-matrix correlators. 
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From technical point of view the road is now open for search of two different generahzations: to 
Dijkgraaf-Vafa models |50il51j . which are not fully specified by the Virasoro constraints alone and rely 
upon intriguing and under-developed theory of check- operators [l9], and to more interesting unitary- 
matrix models with Itzykson-Zuber measures and further to Kazakov-Migdal multi-matrix models 
^4j-[37j, important both for Yang-Mills theory and for the theory of integer partitions. Putting all 
these very different problems into the same context, moreover, underlined by the well established 
theory of free fields on Riemann surfaces [58], is a challenging and a promising perspective. 

Another, but, perhaps, related, open problem is direct derivation of decomposition formula pi.59P 
from integral representations of all the models, bypassing the Virasoro constraints and D-module 
representations. Note that this kind of problem remains unsolved even for the crucially important 
decomposition Zh = U{Zk®Zk)-, describing the double-scaling continuum limit of Hermitian matrix 
model. 
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